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We show that stars in the inner parsec of the Milky Way can be significantly affected by dark matter 
annihilation, producing population-level effects that are visible in a Hertzsprung-Russell (HR) dia- 
gram. We establish the dark HR diagram, where stars lie on a new stable dark main sequence with 
similar luminosities, but lower temperatures, than the standard main sequence. The dark matter 
density in these stars continuously replenishes, granting these stars immortality and solving multiple 
stellar anomalies. Upcoming telescopes could detect the dark main sequence, offering a new dark 


matter discovery avenue. 


“y TARS EVOLVE along well-established tracks within 
k pý which an equilibrium between the gravitational 

Z force and the pressure from nuclear fusion can 
be maintained. Their blackbody spectrum also directly 
links their luminosity, temperature and radius. These re- 
lationships are apparent on the Hertzprung-Russell (HR) 
diagram, which plots the observable stellar luminosities 
and temperatures, tracing various mass-dependent tracks 
stars travel across over their lifetime. 

This standard picture of stellar evolution assumes that 
nuclear fusion is the only process capable of providing 
pressure support for the star. However, in the early 
Universe [l-8], or in regions very close to the Galac- 
tic Center [9-12], dark matter annihilation may pro- 
vide a considerable additional energy source. In ex- 
treme cases, dark matter annihilation can cause stars 
to move backwards along the HR diagrams towards 
proto-stellar configurations (Hayashi track) that have ex- 
tremely large radii and low temperatures compared to 
standard stars [11]. In some cases, dark matter annihila- 
tion can cause stars to be disrupted entirely [13]. 

In this letter we produce, for the first time, a 
population-level description of main sequence evolution 
under the influence of significant dark matter energy in- 
jection — and construct a novel dark main sequence. We 
discover a new branch of higher-mass dark matter pow- 
ered stars that are observationally distinct from the stan- 
dard main sequence. The stability of these stars is anal- 
ogous to the hydrostatic equilibrium of standard stars — 
as the star expands the dark matter capture rate can de- 
crease (see App. |), causing the stars to contract again. 
The continuous accretion of dark matter fuel allows these 
stars to maintain equilibrium eternally, causing the dark 
main sequence to dominate the stellar population. 

Figure | shows our main result: the new dark main 
sequence (top), which characterizes stars very close 
(< 0.1 pc) to the Galactic Center. Compared to the stan- 
dard main sequence (bottom), there are striking differ- 


* isabelle.john@fysik.su.se, ORCID: 0000-0003-2550-7038 
t rleane@slac.stanford.edu, ORCID: 0000-0002-1287-8780 
t linden@fysik.su.se, ORCID: 0000-0001-9888-0971 


Number of Stars 


Dark Main Sequence 
Benchmark 1 


Number of Stars 


Standard Main Sequence 


3x 104 104 6x 103 
Terr [K] 


FIG. 1. The dark main sequence (top) compared to the stan- 
dard main sequence for stellar evolution (bottom), which pre- 
dicts qualitatively different stellar populations in the Galac- 
tic Center. The dark main sequence is shown for benchmark 
parameters of 1 GeV dark matter mass, scattering cross sec- 
tion of 3 x 107°? cm?, and an NFW dark matter profile with 
y = 1.5. 


ences. First, dark matter powered stars can exist in sta- 
ble equilibria at similar luminosities, but lower temper- 
atures (near the Henyey track) than the standard main 
sequence. Second, the dark main sequence lacks low-mass 


stars, as they are overpowered by dark matter and move 
backwards along the Hayashi track [11] — potentially to 
the point that they are entirely disrupted [13]. 


Studying the impact of dark matter on stars at the 
population level allows us to statistically distinguish the 
dark main sequence from the standard main sequence. 
To date, the number of known stars in the inner parsec 
is limited, and more stars are needed to robustly detect 
any effect from dark matter. This will be possible with 
upcoming thirty-meter class telescopes, which will mea- 
sure many new stars near the Galactic Center [14]. 


Interestingly, our initial data regarding stars near the 
Galactic Center already unveils numerous anomalies 
that might be solved by the dark main sequence. Star 
formation models suggest that stars cannot form within 
~ 0.1 pe of the central black hole, where the S-cluster 
stars are found [15]. Rather, the stars must have 
formed elsewhere and migrated towards the Galactic 
Center [16]. Conversely, observations suggest that stars 
in this region are young (S 15 Myr), which indicates 
that the stars might have formed more locally [17]. The 
tension between these two observations is known as the 
“paradox of youth” problem: Galactic Center stars are 
as bright as young stars but show spectroscopic features 
of more evolved stars [16]. Additionally, there is a “co- 
nundrum of old age”, which describes the paucity of old 
and evolved stars, despite the fact that models predict 
the detection of many evolved stars [18, 19]. Finally, 
observations indicate that the initial mass function of 
stars near the Galactic Center is surprisingly top-heavy 
(biased towards high-mass stars), with observations 
indicating that the stellar birth distribution lies between 
M-1*7 [20] and M~°-4° [21], compared to a standard 
Salpeter index of M7235 [22]. These anomalous obser- 
vations are unique to stars at the Galactic Center, and 
are not well explained by conventional stellar evolution 
models. The very inner Galaxy is exactly where the 
dark matter density is high enough for dark matter 
annihilation to substantially replace nuclear fusion as 
the stellar energy source [11, 13, 23-25], allowing stars 
to stay forever young despite their advanced age [26]. 
This explains both the paradox of youth (because 
the stars have had sufficient time to migrate towards 
the Galactic Center), and the conundrum of old age 
(because the stars continue to reside on or near the main 
sequence). Because this effect is maximized in relatively 
high-mass stars (25 Mọ), while less massive stars tend 
to evolve backwards along the Hayashi track or be 
entirely disrupted, the observed stellar distribution will 
naturally have a top-heavy initial mass function [21]. 
These otherwise anomalous observations of Galactic 
Center stars are therefore all naturally predicted on our 
dark main sequence. 


Stellar Population Inputs.— To construct our stan- 
dard and dark HR diagrams, we generate a mock stellar 
population that evolves both with and without the influ- 
ence of dark matter. Throughout this letter, we choose 


relatively standard star formation parameters which are 
also motivated by the observation of S-cluster stars ob- 
served near the Galactic Center [17, 27]. 

We start by randomly selecting stars following a 
Salpeter initial mass function [22] 


dN z 
T ae M; 253 (1) 


where N is the number of stars and M, is the stellar 
mass. We consider stars with masses ranging from 1 to 
20 Mo with a step of 0.05 Mo. 

We set the stellar distances to follow a radial dis- 
tribution that falls as r~? [15, 28, 29], ranging from 
4 x 1075 to 4 x 107? pe. Motivated by the migration sce- 
nario, where stars cannot form closely to the inner black 
hole [16], our stellar evolution models in dark matter sce- 
narios do not include any dark matter energy injection 
until stars enter the main sequence, meaning that the ra- 
dial binning of stars is irrelevant during their pre-main 
sequence life. We finally assume that star formation and 
migration occurs in steady state throughout the 10 Gyr 
lifetime of our simulations, matching the age of the Milky 
Way. 

For simplicity, we assume that the stellar orbits are 
circular, i.e. the dark matter capture rate is constant 
along the orbit. While the orbits of Galactic Center 
stars are highly eccentric, the equilibrium timescale 
for dark matter capture and annihilation is much 
longer than the stars’ orbital period, meaning that 
the effective dark matter annihilation rate is only 
set by the time-averaged orbital radius. Thus, this 
only shifts the exact capture rate for each star that is 
found somewhere along the region around Sagittarius A* 
that we consider, and does not largely change our results. 


Dark Matter Inputs.— We generate two dark matter 
model benchmarks. Qualitatively, we find similar behav- 
ior across both, and show that similar behavior is ex- 
pected for a wide range of dark matter parameters. In 
both cases, we assume a dark matter mass of 1 GeV, but 
vary the following parameters: 

Benchmark 1. We set the dark matter-nucleon scat- 
tering cross section to oyy = 3 x 1075 cm?. For this, 
we follow our constraints on the scattering cross section 
from Ref. [13] for the star S2, which stands as the best- 
measured S-cluster object. To set the dark matter den- 
sity along the different orbits, we choose a generalized 
NFW profile [30] with index y = 1.5. This corresponds 
to dark matter densities of ~ 10° — 101° GeV/cm? for 
the distances considered here. 

Benchmark 2. We choose a smaller dark matter- 
nucleon scattering cross section of oyy = 1 x 107%" cm?. 
For the dark matter density, we choose a standard NF W 
(y = 1.0), but include a density spike that arises from the 
adiabatic accretion of matter by Sagittarius A* [31, 32]. 
We follow the spike model from [13, 33], i.e. y = 1.0, 
Rspike = 10 pe and Yepike = 7/3. This results in higher 
dark matter densities compared to Benchmark 1, and 
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FIG. 2. The impact of dark matter burning on stars over a range of masses and Galactic Center distances, for Benchmark 
Models 1 (left panel) and 2 (right panel). Green shades indicate how much main sequence stellar evolution is slowed down by 
the added dark matter power compared to the standard evolution scenario. This corresponds to stars that are more massive and 
experience a lower dark matter density. Stars in higher density regions are significantly changed and show atypical evolution: 
The stars maintain equilibrium through dark matter burning and only rely on small amounts of fusion, making them reside 
on the Henyey track (yellow shade), or have no fusion at all, which corresponds to the Hayashi track (orange shades). This 
essentially “freezes” their evolution, as stellar evolution is driven by nuclear fusion. The continuous re-supply of captured dark 
matter acts as an effectively infinite power source, making these stars immortal. 


the density increases more steeply towards the Galactic 
Center, i.e., ~ 101? — 1016 GeV/cm?. 

Our benchmarks are not excluded by direct detection 
for spin-dependent dark matter models [34-38], and 
even given improved direct detection limits in the future, 
there are a wide range of additional parameter points 
that provide qualitatively similar results to those we 
show here. For example, our qualitative results are not 
affected for lower dark matter masses. In addition, at 
lower cross sections, large dark matter densities can 
compensate for any decrease in the signal. Finally, 
we note that for extremely high densities, the efficient 
annihilation of dark matter may cause the dark matter 
density to saturate at a value that depends on the dark 
matter mass and annihilation cross section. However, in 
the case of models that use stellar capture to focus dark 
matter annihilation, our scattering cross-section limits 
continue to apply even for extremely small annihilation 
cross-sections that will not saturate up to very large 
densities. 


Stellar Evolution Simulations. — We use the stellar 
evolution code MESA [39-44], version r24.03.1 with 
Software Development Kit (SDK) version 23.7.3, 
to investigate the effect of dark matter energy injec- 
tion on a wide range of stellar masses. We base our 
model on the default work model, but implement a 
custom module to calculate the capture rate at each 
time step according to the current stellar radius. To 
obtain the capture rate, we follow the formalism of 
Ref. [45], under the assumption of capture-annihilation 


equilibrium, and add the annihilation energy into the 
stellar core through the inject_uniform_extra_heat 
parameter. The injection of the additional energy can 
be restricted to a certain mass region of the star through 
the parameters min-q-for-uniform-extra-heat and 
max_q_for_uniform_extra_heat. We consider the stellar 
core to correspond to the inner 10% of the star, but 
have tested that injecting the dark matter annihilation 
energy in the inner 1% only results in slight numerical 
differences, while uniform injection throughout the star 
only changes our results by a factor of a few, which does 
not have qualitative effects. We assume that all stars 
initially form and evolve in a low dark matter density 
environment until their early main sequence phase, 
i.e. their pre-main sequence evolution is not affected 
by dark matter. We assume that these stars migrate 
as main sequence stars towards the Galactic Center 
where their subsequent evolution is influenced by dark 
matter annihilation. As discussed above, this migration 
scenario is a commonly considered case, due to the fact 
that stellar formation close to Sagittarius A* is expected 
to be inhibited [16]. In MESA, we start injecting the 
additional energy when the hydrogen fraction in the star 
center drops below 0.67, which corresponds to an early 
main sequence core composition, through the parameter 
xa_central_lower_limit. Once the star experiences 
dark matter energy injection, we evolve the star until it 
reaches the red giant branch phase, using the stopping 
condition stop_at_phase_ZACHeB, or until an age of 
10 Gyr is reached, corresponding to the age of the Galaxy. 


Results and Discussion.— Figure 2 shows the ef- 
fects of dark matter energy injection on stars of differ- 
ent masses and distances from the Galactic Center for 
Benchmark 1 and 2. These models only differ in the 
dark matter-nucleon scattering cross section and the un- 
derlying dark matter profile, that determines the amount 
of dark matter a star can capture along its orbit. We find 
that stellar evolution is affected in various ways: 


(1) Longer-lived. The green shaded regions correspond 
to stars that are not significantly or only moderately 
changed by dark matter burning. Their evolution 
follows the typical stages, but as part of the nuclear 
fusion power needed to maintain hydrostatic equilib- 
rium is replaced by dark matter power, hydrogen in 
the core fuses more slowly and thus the stars’ evo- 
lution is slowed down. This makes the stars longer- 
lived. While light green shades correspond to stars 
that are only slightly affected as their evolution times 
are similar to the standard scenario, dark green re- 
gions represent stars that can become as much as 
10-100 longer-lived. 


(2) Stuck on Henyey track. The yellow region corre- 
sponds to stars that are found on the Henyey track 
in the HR diagram, which, in the standard evolu- 
tionary scenario, is a rapid phase just before stars 
enter the main sequence, when hydrogen burning is 
ignited and increasing. This means that at the onset 
of dark matter burning in the stellar core, stars can 
move “backwards” in evolution, returning from the 
main sequence to a stage with only small amounts 
of nuclear fusion, and maintain equilibrium primar- 
ily through dark matter burning. Dark stellar evo- 
lution on the Henyey track is extremely slow, until 
eventually the star does run out of hydrogen in the 
core and moves back to a more typical evolution on 
the red giant track. Thus, these stars are in princi- 
ple extremely long-lived and would eventually follow 
the standard evolution stages, but, as in our simula- 
tions we do not evolve stars for longer than 10 Gyr, 
which corresponds to the Galaxy’s age, these stars 
never move off the Henyey track. Note that many of 
the stars in scenario (1) that are significantly longer- 
lived (Z 5 times) also spend some time on the Henyey 
track, but experience a lower amount of dark matter 
burning and thus need to undergo more nuclear fu- 
sion, forcing them to move back to the main sequence 
and resume evolution earlier. 


(3) Stuck on Hayashi track. The orange region corre- 
sponds to dark matter powered stars that are found 
in the region of the Hayashi track in the HR diagram, 
which, in the standard picture, describes pre-main se- 
quence stars that are still contracting and have not 
yet started nuclear fusion. At the onset of dark mat- 
ter burning, main sequence stars evolve “backwards”, 
similar to the stars halting on the Henyey track de- 
scribed above. In this scenario however, dark mat- 
ter burning replaces nuclear fusion entirely, and stars 
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FIG. 3. Dark HR diagram for a population of Galactic Center 
stars, for our dark matter benchmark 2, with 1-GeV dark 
matter mass, scattering cross section of 107°" cm?, and NFW 
dark matter profile with a spike (see text for details). 


move all the way to the Hayashi track. The stars do 
not undergo any further evolutionary stages as long 
as the dark matter capture rate remains sufficiently 
high to prevent nuclear fusion from resuming. A star 
like this might look like a young, still-forming star 
but has features of a star that has undergone nuclear 
fusion in the past — and is effectively immortal. 


(4) Beyond Hayashi track. The dark orange region corre- 
sponds to stars moving back onto the Hayashi track, 
similar to case (3). However, in this case, the stars 
move further up the Hayashi track to even earlier 
evolutionary states. In our MESA simulations, these 
stars end up higher on the Hayashi track than their 
starting position. We emphasize that this scenario 
is highly reliant on the modelling methods (as MESA 
does not simulate of the very early stages of star for- 
mation). It is possible that these stars could expand 
even further or be entirely disrupted by the dark mat- 
ter energy injection. A scenario like this may be con- 
sistent with G-objects, a recently discovered strange 
class of objects found at similar orbits around Sgr 
A* as the S-cluster stars that resemble clouds of gas 
surrounding a denser core [46, 47]. 


In both benchmarks, massive stars further out from the 
Galactic Center are the least affected, while light stars 
can become immortal even at distances where the dark 
matter density is not sufficient to significantly affect more 
massive stars. The different dark matter density profiles 
are also reflected in Fig. 2: While Benchmark 1 (left 
panel) uses an NFW profile, Benchmark 2 (right panel) 
includes a dark matter spike that steepens the profile 
and, accordingly, the region of the parameter space in 
which stars are immortal. 

Figure | shows HR diagrams for a population of Galac- 
tic Center stars with and without dark matter for Bench- 
mark 1, and Figure 3 shows the dark HR diagram for the 


our additional Benchmark 2. We include evolutionary 
stages from the early main sequence up until the stars 
leave the main sequence and start hydrogen shell burn- 
ing. We do not include stars that have become red giants 
and start core Helium burning, i.e. most of the stars on 
the HR diagram are main sequence stars. Most stars 
lie on the main sequence, i.e. the diagonal line, where 
more massive stars are at the top left, as they are hotter 
and more luminous, and lighter stars are at the bottom 
right. The small number of stars scattered to the right 
of the main sequence are stars that are older stars that 
have just moved away from the main sequence towards 
the red giant branch; which manifests in a lower temper- 
ature compared to the main sequence. 


For our dark main sequence in Figs. | and 3, the ef- 
fect of dark matter burning results in strikingly distinct 
features. Compared to the standard main sequence, the 
Hayashi tracks (approximately vertical tracks at low tem- 
peratures) are densely populated with immortal stars 
(see scenarios (3) and (4) above), and the Henyey tracks 
(approximately horizontal tracks at high luminosities), 
are clearly visible (see scenario (2) above), containing 
stars that are mostly powered by dark matter burning 
with only small amounts of nuclear fusion, making their 
evolution extremely slow. Additionally, the region where 
light main sequence stars are located (lower luminosi- 
ties and temperatures), are only sparsely populated, as 
most of these stars have moved to the Hayashi tracks. 
In Benchmark 2 (Fig. 3), more of these low-mass stars 
remain on the main sequence compared to Benchmark 1 
(Fig. 1), due to the steeply increasing dark matter density 
of the spike model, as can also be seen in Fig. 2. We note 
that the striped features observed in the Henyey branch 
are due to our distance (and thus dark matter density) 
binning, which sets the stellar luminosity, and are thus 
not a physical feature. 


Finally, we note that the fraction of high-mass stars 
increases in dark matter scenarios due to the effect 
of stellar immortality. If only standard stellar models 
were considered, the large fraction of high-mass stars 
would be misconstrued as evidence of a top-heavy initial 
mass function (IMF). Intriguingly, observations have 
found evidence of a top-heavy IMFs at the Galactic 
center [20, 21, 48]. We find that dark matter would 
harden the IMF, compared to a Salpeter IMF, by 
approximately ~ 0.5 — 2.5, as shown in Figs. | and 3. 


Summary and Outlook.— We have produced the first 
population-level prediction for stellar evolution under the 
presence of dark matter capture and annihilation near the 
Galactic Center. We have established a new “dark main 
sequence” that is observable in the Hertzsprung-Russell 
diagram of stars observed in dark matter rich environ- 
ments. This sequence has three observable characteristics 
that clearly separate it from standard main sequence stel- 
lar evolution: (1) a top-heavy initial mass function along 
the main sequence branch, produced by the efficient re- 
moval of low-mass stars off of the main sequence and into 


locations on the Hayashi or Beyond-Hayashi tracks, (2) 
a bright and densely populated Hayashi track, produced 
by low-mass main sequence stars which achieve eternally 
stable equilibria in configurations with large stellar radii 
and low stellar temperatures, (3) a bright and densely 
populated Henyey track at higher-masses, populated by 
effectively immortal (typical lifespans >10 Gyr) high- 
mass stars with similar luminosities but slightly lower 
temperatures than their main sequence configurations. 

Excitingly, there are several anomalies in Galactic Cen- 
ter stellar observations which can be simultaneously ex- 
plained by stellar evolution along our dark main se- 
quence. The top-heavy initial mass function near the 
Galactic Center is automatically produced by our sim- 
ulations with dark matter, while both the paradox of 
youth and conundrum of old age can be explained if rel- 
atively high-mass stars are powered indefinitely by dark 
matter annihilation. Finally, observations of G-objects 
(extremely large gas clouds which may hold a central 
star), may be consistent with Beyond-Hayashi evolution, 
although the stability of such solutions falls beyond the 
parameters of our current MESA modeling. We will inves- 
tigate this in upcoming work. 

Going forward, our stellar population predictions 
have important implications for the extremely difficult 
observations of Galactic Center stars. The vastly cooler 
temperatures of Hayashi track stars may impede their 
detection unless observational efforts specifically look for 
such objects. Additionally, the dark Henyey track stars 
may be mischaracterized as standard main sequence 
stars if the possibility of Henyey track evolution is not 
considered. Verifying or ruling out dark main sequence 
evolution will require the detection of new stars for 
sufficient sampling and to sufficiently populate the 
dark main sequence. This will be possible with new 
telescopes, such as thirty-meter class telescopes, or with 
additional stars measured with the Very Large Telescope 
(VLT) or the Keck Observatory. 
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I. DARK MATTER CAPTURE RATES FOR ACTIVELY EVOLVING STARS 
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FIG. S1. Left panel: Dark matter capture rate as a function of stellar radius, for a scattering cross section of 10737 cm?. Right 
panel: Capture rate as a function of scattering cross section, for stellar radii of 1 and 10 solar radii, and a fixed stellar mass of 
1 solar mass. The dashed line shows the transition cross sections for the two stars. To calculate these capture rates, we assume 
an example dark matter mass of 0.1 GeV, orbital velocity of 10 cm/s and dark matter density 2 x 10’? GeV/cm?. 


Figure 51 (left panel) shows the dark matter capture rates expected as the star undergoes evolution, i.e. taking into 
account the changing mass and radius of the star. The first point to note, as highlighted by the left figure, is that 
the capture rate changes as the stellar radius changes, and so it is important to take the updated stellar radius into 
account for the capture rate as the star evolves. This implementation is an improvement over previous studies [13]. 
The second point is that for a fixed stellar mass and fixed dark matter scattering cross section, the dark matter 
capture rate can decrease as the stellar radius increases. The reason for this can be understood from the right figure. 

Figure 51 (right panel) shows the capture rate as a function of scattering cross section, for two example stars of 
radius 1 and 10 times the solar radius, for a fixed stellar mass of 1 solar mass. We also show as the dashed line the 
transition cross section for these two example stars. Here, the transition cross section is Ctr = 7R?/N, where N is 
the number of stellar nucleons available for scattering. Physically, this corresponds to a mean free path of the dark 
matter of about the size of the star, and therefore corresponds to the cross section above which multiple scatters 
can be needed for capture to occur. Below this cross section, dark matter can only be captured with roughly one 
scatter, as the scattering is so weak that on average one or zero scatters occur. As the mass of these two example 
stars stays the same, the number of targets also remains the same, except they are distributed over a larger volume, 
the transition cross section between the two simply scales as their radius ratio squared. 

As Fig. 51 (right) shows, larger radii stars do capture more dark matter in the sufficiently-high cross section regime 
where roughly all incoming dark matter is captured, as their surface areas available for capture are larger. However, 
less dense objects have lower escape velocities, so require more scatters (i.e. a larger cross section) to bring the dark 
matter to a low enough speed to be captured, which is why the larger star needs a larger cross section increase above 
its transition cross section at ~ 10~3° cm? to reach its geometric capture rate. In comparison, the smaller radius 
star increases less above its ~ 10735 cm? transition before reaching geometric capture, because it is more efficient at 
capture (and so needs fewer scatters) due to its higher escape velocity. Once comparing the capture rates for these two 


objects at their transition cross section, the puffier star is already down a larger penalty compared to its geometric 
maximum capture than the more dense star, which is then cancelled out by its size, such that at their respective 
transition cross sections (which are two orders of magnitude apart) they actually have the same capture rate. This 
means that, because the capture rate decreases linearly with cross section below the transition cross section, the 
puffier, larger radius star has a lower capture rate at smaller cross sections, as its transition cross section is higher to 
start with. This can lead to interesting effects during the stellar evolution process as the star becomes larger, as the 
dark matter capture and annihilation driving the growth in radii consequently decreases. 
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